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ABSTRACT 


The  Naval  Postgraduate  School  open  circuit  oscillating 
flow  wind  tunnel  was  used  to  study  the  blowing  requirements 
to  maintain  an  attached  turbulent  boundary  layer  in  an  oscil- 
lating freestream  flow  with  an  adverse  pressure  gradient. 
Boundary  layer  separation  was  visualized  through  the  use  of 
tufts.   Freestream  flow  oscillation  frequency  was  found  to 
have  an  effect  on  the  flowing  required  to  maintain  attached 
flow.   This  frequency  dependence  exhibits  characteristics 
which  suggest  resonant  behavior. 
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I.   INTRODUCTION 

The  study  of  viscous  unsteady  flows  have  become  increas- 
ingly important  in  recent  years.   Typical  examples  of  un- 
steady flows  are  found  in  helicopter  blade  aerodynamics  and 
in  gas  turbine  engine  internal  flows.   The  development  of 
advanced  design  concepts  in  helicopters,  ships,  and  propul- 
sion systems  requires  a  basic  understanding  of  the  nature  of 
unsteady  flows.   The  complexity  of  the  problem  does  not  lend 
itself  readily  to  either  analytic  or  experimental  study, 
even  when  limited  to  flow  with  laminar  boundary  layers . 
Furthermore,  most  of  the  problems  of  interest  are  characterized 
by  turbulence  and  turbulent  boundary  layers .  which  greatly 
increases  the  difficulty  of  the  study.   Regardless  of  the 
difficulties  involved  the  field  is  of  such  importance  as  to 
warrant  the  expenditure  of  time  and  effort  to  widen  the  un- 
derstanding of  its  characteristics. 

The  present  study  was  to  experimentally  examine  the  ef- 
fect of  different  levels  of  blowing  on  turbulent  boundary 
layer  separation  in  an  oscillating  freestream.   The  effects 
of  changes  in  frequency  and  magnitude  of  the  freestream 
oscillation  on  the  blowing  requirements  to  maintain  an  at- 
tached boundary  layer  were  studied. 
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II.   BACKGROUND 

A.   NON-STEADY  FLOWS 

1 .   Analytic  Studies 

Previous  analytic  studies  have  been  confined  mostly 
to  flows  which  have  a  laminar  boundary  layer.   The  earliest 
was  Stokes  [Ref.  1]  study  of  the  doubly  infinite  flat  plate 
oscillating  in  its  own  plane  in  a  fluid  at  rest.   Rayleigh 
[Ref.  2]  considered  the  second  order  effects  of  this  prob- 
lem.  Schlichting  [Ref.  3]  expanded  Rayleigh's  work  to  in- 
clude the  boundary  layer  assumptions.   All  of  these  were 
simple  unsteady  flows  without  a  mean  flow  or  pressure  gra- 
dient.  Lighthill  [Ref.  4]  later  treated  the  case  of  small 
magnitude,  low  frequency,  sinusoidal  oscillations  super- 
imposed on  a  steady  flow  about  a  cylinder  and  a  semi-infinite 
flat  plate.   Using  small  perturbation  theory  and  restricting 
the  solution  to  first  order  terms  he  determined  that ,  at  low 
frequencies  the  flow  is  essentially  quasi-steady.   The  un- 
steady flow  had  the  same  characteristics  for  any  instantan- 
eous magnitude  of  freestream  velocity  as  a  steady  flow  of 
the  same  velocity.   As  frequency  is  increased  a  limiting 
frequency  is  reached  beyond  which  the  boundary  layer  reacts 
as  it  would  without  a  mean  flow.   Lin  [Ref.  5]  analyzed  the 
high  frequency  case  for  large  amplitude  oscillations.   He 
found  that  the  governing  equations  became  essentially  lin- 
ear at  high  frequencies  allowinp  separation  of  oscillatory 
and  steady  components.   Nickerson  [Ref.  6]  expanded  this 
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work  to  include  higher  order  perturbation  terms.   Barriol 
and  Lucius  [Ref.  7]  used  numerical  methods  to  obtain 
asymptotic  solutions  to  the  boundary  layer  equations  for 
oscillating  flow  on  a  semi-infinite  flat  plate  with  no 
pressure  gradient.   Their  solution  agrees  with  those  of 
Lighthill  and  Nicker son. 

The  previously  mentioned  studies  have  all  dealt  with 
laminar  boundary  layer  flows .   The  analytic  approach  to  tur- 
bulent boundary  layer  flow  has  apparently  been  less  fruitful. 
Karlssen  [Ref.  8]  separated  the  fluctuations  into  periodic 
and  random  components ,  and  by  averaging  over  a  complete 
period  obtained  equations  similar  to  the  steady  turbulent 
boundary  layer  equations,  but  with  an  additional  fluctuating 
term.   Other  analytic  works  [Ref.  9,  10]  have  used  linearized 
solutions  to  solve  the  pressure  distribution  on  an  airfoil 
at  low  angle  of  attack,  but  are  applicable  to  only  a  small 
number  of  cases. 

2 .   Experimental  Studies 

A  great  deal  more  work  has  been  done  experimentally 
than  analytically,  but  there  are  still  large  gaps  in  the 
field.   Nickerson  [Ref.  6]  partially  verified  both  his 
analysis  and  that  of  Lighthill.   Miller  [Ref.  11]  studied  ■ 
the  transition  phenomena  on  a  flat  plate  in  oscillating  flow. 
He  was  able  to  confirm  some  of  the  previous  analytic  pre- 
dictions.  In  addition  he  determined  the  transition  Reynolds 
number,  turbulent  Reynolds  number,  and  turbulent  intermit- 
tency  factor  for  oscillating  Blasius  type  flow.   Despard 
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[Ref .  12]  investigated  the  separation  of  a  laminar  boundary 
layer  in  oscillating  flow.   He  proposed  that  the  definition 
of  separation  for  a  laminar  boundary  layer  in  an  oscillating 
freestream  be  the  initial  occurrence  of  zero  velocity  or 
reverse  flow  at  some  point  in  the  velocity  profile  through- 
out the  entire  cycle  of  oscillation.   He  was  also  able  to 
make  some  prediction  about  the  behavior  of  the  separation 
point . 

Morrissey  [Ref.  13]  investigated  the  effects  of 
large  amplitude  flow  oscillations  on  the  heat  transfer  from 
a  flat  plate  with  a  turbulent  boundary  layer.   Jacobs  [Ref. 
14]  studied  the  effect  of  oscillating  mean  freestream  on 
the  turbulent  intensity  distribution  in  a  turbulent  boundary 
layer.   These  two  studies  indicate  that  unsteady  flows  ex- 
hibit no  significant  alteration  in  the  character  of  the  eddy 
diffusivity  or  the  turbulent  intensity  distribution  when 
compared  to  a  steady  turbulent  flow.   Banning  [Ref.  15]  in- 
vestigated the  pressure  distribution  on  an  airfoil  in  a 
turbulent  oscillating  freestream.   Others  have  studied  the 
effects  of  oscillating  flows  with  compressibility  effects 
included . 

B.   BLOWING  TO  AVOID  SEPARATION 
1 .   Experimental  Studies 

Boundary  layer  control  through  blowing  has  been  with 
us  almost  as  long  as  the  boundary  layer  concept  itself.  The 
idea  arose  initially  from  the  use  of  slots  as  a  lift  augmen- 
tation device  ,  which  was  suggested  by  Lachmann  of  Germany  in 
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1918  and  later  tried  by  him  and  Handley-Page  in  England. 
Betz  [Ref.  16]  theorized  that  the  effect  of  the  slots  was 
to  accelerate  the  boundary  layer.   Baumann  [Ref.  17]  was 
led  by  this  interpretation  to  replace  the  air  passing  through 
the  slot  with  air  ejected  from  the  wing  interior.   This  pro- 
duced roughly  the  same  effect  as  the  slots  with  the  added 
advantage  of  being  able  to  control  the  efrect  by  control  of 
the  pressure  inside  the  wing. 

Until  the  1940' s  little  more  was  done  with  blowing 
except  some  experiments  which  used  it  as  an  alternative  for 
ailerons.   This  was  done  by  blowing  out  a  slot  over  a  short 
wing  section  to  induce  high  lift.   By  differential  blowing 
on  the  wing  a  rolling  moment  was  produced  similar  to  that 
of  an  aileron.   In  1942-4  3,  triggered,  no  doubt,  by  the  war, 
experiments  began  to  be  conducted  in  France,  Germany  and  the 
United  States,  followed  shortly  by  Britain.   Probably  the 
most  important  contribution  of  this  era  was  made  in  1948  by 
Poisson-Quinton  [Ref.  18].   Prior  to  this  time  the  common 
parameter  to  measure  blowing  was 

CQ  =  D^S  (II-1) 

x       oo 

where  Q.  is  the  volumetric  flow  through  the  slot,  Uro  the 
mean  freestream  velocity  and  S  the  characteristic  area  of 
interest.   Schwier  [Ref.  19]  had  shown  that  narrow  slots 
were  more  efficient  than  wide  ones.   Poisson-Quinton  verified 
this  and  suggested  that  a  more  appropriate  blowing  parameter 
was 
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m.V. 

where  m .  is  the  mass  flow  of  blown  air,  V.  the  net  velocity 
]  3 

at  the  exit  assuming  isentropic  expansion  from  the  total 
pressure  and  temperature  inside  the  jet  supply  duct,  and  q 
the  freestream  dynamic  pressure.   Subsequent  experimental 
work  supported  the  use  of  this  parameter  ^o  predict  the 
effectiveness  of  blowing. 

Unfortunately  all  the  early  efforts  in  boundary 
layer  control  through  blowing  were  hampered  by  the  great 
difficulties  involved  in  practical  applications.   Major 
obstacles  were  a  supply  of  high  pressure  air  and  the  weight 
penalty  associated  with  the  ducting.   However,  with  the 
advent  of  gas  turbine  engines  for  aircraft  propulsion  a 
good  supply  of  air  became  available.   With  this  development 
and  advances  in  metal  alloys  and  other  lightweight  materials 
of  high  strength,  blowing  has  become  much  more  attractive. 
Since  1950  the  use  of  blowing  has  proceeded  at  an  accelerat- 
ing pace.   Unfortunately  little  work  has  been  accomplished 
in  the  area  with  which  this  study  is  concerned.   The  majority 
of  the  work  has  been  in  the  use  of  blowing  for  lift  augmenta- 
tion by  delaying  stall  or  by  providing  circulation  control.. 
All  of  these  experiments  present  results  in  terms  of  Cp  jCj,, 

and  C   and  their  variation  with  C  .   Also  the  values  of  C 
m  y  y 

required  for  lift  augmentation  are  much  greater  than  that 
required  to  keep  the  boundary  layer  from  separating  under 
the  influence  of  an  adverse  pressure  gradient.   As  a  resu] I 
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there  has  been  little  interest  in  studying  the  amount  of 
blowing  required  to  maintain  streamline  flow.   The  only  in- 
formation to  be  gleaned  from  this  work  was  that  higher  pres- 
sure gradients  require  higher  C   to  overcome  separation. 
2 .   Analytic  Studies 

Theoretically,  most  of  the  interest  lies  in  increas- 
ing lift ,  rather  than  detailed  study  of  blowing  requirements 
to  maintain  attachment.   There  are,  however,  a  few  notable 
exceptions.   Carriere  and  Eichelbrenner  [Ref.  2  3]  developed 
a  method  for  calculating  the  conditions  for  flow  re-attach- 
ment by  a  jet  discharging  against  adverse  pressure  gradient. 
Their  analysis,  unfortunately,  was  heavily  dependent  on  the 
availability  of  empirical  profiles  from  experimentation  and 
several  arbitrary  assumptions.   Within  these  limits  it  did; 
however,  provide  solution  to  cases  with  weak  pressure  gra- 
dients.  Kozlous  and  Zyanyak  [Ref.  24]  have  recently  anal- 
yzed a  laminar  boundary  layer  in  an  unsteady  incoming  motion 
around  a  body  of  arbitrary  shape  with  either  suction  or 
blowing.   They  used  a  six  degree  polynomial  to  describe  the 
boundary  layer  under  these  conditions.   Their  result  was  an 
integral  equation  which  could  be  integrated  directly,  or 
reduced  to  a  quadrature,  depending  on  the  velocity  profile. 
They   were  able  to  verify  their  analysis  for  a  symmetric 
wing  in  a  start-up  situation.   This  can  be  extended  to  an 
oscillating  flow  by  varying  the  description  of  the  incoming 
motion.   Although  derived  for  lift  and  drag  predictions  in 
laminar  flows  with  blowing,  this  analysis  shows  great  promise 


17 


for  being  extended  to  turbulent  boundary  layer  blowing  re- 
quirement predictions.   The  greatest  difficulty  with  this 
extension,  or  any  prediction  of  blowing  in  turbulent 
regimes ,  is  that  much  of  the  energy  imparted  to  the  boun- 
dary layer  is  by  turbulent  mixing  which  defies  accurate 
modeling . 

C.   APPLICATIONS 

Recently  there  has  been  work  accomplished  toward  the 
application  of  boundary  layer  control  in  an  oscillating 
flow.   Englar  and  Williams  [Ref.  12,  23,  24]  have  applied 
boundary  layer  energization  through  blowing  to  augment  the 
lift  of  a  submarine  stern  plane  and  a  symmetric  airfoil  at 
angle  of  attack.   They  have  also  applied  tangential  blowing 
to  the  blades  of  a  helicopter  rotor  system  (the  circulation 
control  rotor)  with  encouraging  results. 

There  has  not,  however,  been  any  detailed  studies  of 
the  actual  flow  over  the  surface  and  its  behavior  under  the 
conditions  of  high  angle  of  attack,  large  pressure  gradient, 
turbulence,  oscillating  freestream,  or  blowing  to  energize 
the  boundary  layer.   The  prediction  of  blowing  requirements 
under  these  conditions  will  for  the  foreseeable  future , 
rely  heavily  on  empirical  methods  based  upon  experimental 
results.   The  present  study  is  an  attempt  to  provide  some 
of  this  experimental  data. 
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III.   DESCRIPTION  OF  APPARATUS 

A.   OSCILLATING  FLOW  WIND  TUNNEL 
1 .   General  Description 

The  experimental  work  was  conducted  in  the  lowspeed, 
oscillating  flow  wind  tunnel  located  in  the  Aeronautics 
Laboratories  of  the  Naval  Postgraduate  School.   This  wind 
tunnel  is  of  open  circuit  design,  with  a  24-inch  square  by 
223-inch  long  test  section.   A  plan  view  of  the  tunnel  is 
presented  in  Figure  III.l.   The  tunnel  inlet  is  eight  feet 
square,  resulting  in  a  16:1  contraction  ratio.   Three  high 
solidity  screens  located  in  the  inlet  section  just  upstream 
of  the  nozzle  produce  measured  freestream  turbulence  in- 
tensities of  0.261  to  0.413  percent  for  the  velocities 
encountered  in  the  present  work. 

The  wind  tunnel  drive  consists  of  two  Joy  Axivane 
Fans  in  series,  each  of  which  has  an  internal,  100  horse- 
power, direct  connected,  17  5  0  rpm  motor.   The  fan  blades 
are  internally  adjustable  through  a  pitch  range  of  25  to 
55  degrees,  providing  a  wide  operating  range  of  test  section 
velocities .   Two  sets  of  variable  inlet  vanes  ,  located  im- 
mediately upstream  of  each  fan,  are  externally  operated  to 
provide  control  of  test  section  velocity.   These  vanes  are 
of  multileaf  design,  and  preswirl  the  air  in  the  direction 
of  fan  rotation  to  reduce  fan  capacity.   The  total  range  of 
tunnel  velocity  is  from  10  to  250  ft.  per  second. 
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2 .   Rotating  Shutter  Valve 

Two  fundamental  methods  of  creating  an  oscillating 
flow  environment  have  been  employed  in  the  past.   Nickerson 
[Ref.  6]  introduced  oscillations  by  oscillating  the  model 
in  a  steady  flow  environment.   This  method  severely  re- 
stricts the  range  of  attainable  frequencies  because  of 
mechanical  complications ,  and  also  introduces  measurement 
difficulties.   The  other  method  is  to  actually  oscillate  the 
flow  over  a  stationary  model.   Hill  [Ref.  25]  used  a  sliding 
shutter  to  impose  oscillations  on  the  freestream  but  was  re- 
stricted by  mechanical  limitations  to  low  frequencies. 

The  most  successful  method  of  obtaining  an  oscillat- 
ing flow  with  large  ranges  of  frequency  and  amplitude  was 
that  employed  by  Karlssen  [Ref.  8]  and  later  by  Miller  [Ref. 
11]  in  his  investigation  of  transition  phenomena.   A  rotating 
shutter  valve,  immediately  downstream  of  the  test  section, 
is  employed  to  superimpose  a  periodic  variation  of  velocity 
on  the  mean  flow.   The  method  used  in  the  present  investiga- 
tion is  identical  to  that  employed  by  Miller.   The  shutter 
valve  consists  of  four  horizontal  steel  shafts  equally  spaced 
across  the  test  section.   The  shafts  are  slotted  to  accommo- 
date flat  blades  of  various  widths  ,  forming  a  set  of  four  ■ 
butterfly  valves  spanning  the  test  section.   Figure  III. 2  is 
a  photograph  of  the  shutter  valve .   Each  blade  is  driven  from 
its  immediate  neighbor  by  means  of  a  timing  belt  and  pulley 
arrangement.   The  bottom  shaft  is  driven  by  a  five  horse- 
power, variable  speed,  electric  motor,  through  a  timing  belt 
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and  pulley.   An  intermediate  shaft  between  the  motor  and 
shutter  valve  permits  a  wide  variety  of  pulley  ratios.   This 
drive  arrangement  provides  a  frequency  range  from  0 . 4  cycles 
per  second  to  the  first  critical  frequency  of  933  cycles 
per  second.   The  electric  motor  presently  in  use,  however, 
restricts  the  oscillation  frequency  to  a  maximum  of  240 
cycles  per  second.   The  amplitude  of  oscillation  is  controlled 
by  blade  width.   Test  section  closure  may  be  varied  from  25 
to  100  percent.   The  resulting  amplitude  of  oscillation  of 
test  section  velocity  is  a  function  of  frequency,  mean 
velocity,  and  pressure  gradient.   In  this  investigation 
blades  producing  83.3  percent  closure  were  used  resulting  in 
a  perturbation  range  from  5  to  2  5  percent  of  the  local  mean 
freestream  velocity.   Figure  III . 3  is  a  picture  of  the  shut- 
ter valve  drive  arrangement. 
3 .   Test  Section 

The  wind  tunnel  test  section  is  shown  in  Figure  111.4. 
Continuous  pieces  of  two-inch  thick  aluminum,  24  inches  wide 
and  22  3  inches  long,  form  the  upper  and  lower   test  section 
walls.   Each  of  the  sidewalls  consists  of  three  two-inch 
thick  panels  of  stress  relieved  lucite.   For  this  investiga- 
tion the  central  sidewall  panel  on  the  opposite  side  of  the 
tunnel  from  the  control  console  was  replaced  with  two-inch 
thick  plywood  to  facilitate  the  mounting  of  instrumentation 
and  access  to  the  model  plenum.   The  Lucite  panels  on  the 
console  side  of  the  test  section  are  hinged  and  may  be 
raised  hydraulically ,  providing  access  to  the  test  section. 
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Figure  III. 5  shows  the  test  section  with  the  door  open. 
The  heavy  construction  of  the  test  section  is  intended  to 
minimize  deflections  induced  by  rapid  changes  in  static 
pressure . 

Previous  test  on  the  test  section  velocity  profile 
have  shown  the  velocity  variation  is  less  than  one  percent 
of  mean  to  within  three  inches  of  the  wall  [Ref.  12]. 
Figure  III.  6  is  a  photographic  view  of  the  test  section  and 
control  console  with  the  model  in  place. 

B.   MODEL 

The  model  used  in  this  investigation  consisted  of  two 
plugs  placed  in  the  test  section  with  a  flat  plate  halfway 
between  them.   Each  plug  had  a  smoothly  curved  leading  edge, 
followed  by  42-inch  straight  section,  then  a  diverging  sec- 
tion as  shown  in  Fig.  III.  7.   The  maximum  thickness  of  each 
plug  was  six  inches  and  the  plate  was  \   inch  leaving  5-3/4 
inches  between  each  plug  and  the  plate.   The  diverging  sec- 
tion departs  from  the  horizontal  by  27°.   The  slant  length 
was  13  inches  and  the  characteristic  area  chosen  for  calcu- 
lation of  blowing  coefficient  was  the  plane  area  of  the 
diverging  section  which  was  305.5  square  inches. 

Each  plug  contained  a  1700  cubic  inch  plenum  chamber. 
The  entrance  to  the  plenum  was  a  two-inch  diameter1  hole  to 
which  the  blowing  air  supply  hose  was  connected.   The  blow- 
ing slots  were  2  3.5  inches  wide .   The  upper  slot  height  was 
.041  inch  and  the  lower  was  .055  inch.   The  difference  in 
slot  heights  provided  information  on  the  blowing  coef f ic \ 
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requirements  as  a  function  of  slot  height.   Because  of  its 
smaller  slot  height  the  upper  jet  wo.uld  need  less  mass  flow 
but  higher  jet  velocity  to  achieve  the  same  blowing  coeffi- 
cient as  the  lower  slot.   Figure  III. 8  shows  the  plenum  and 
slot  configuration.   Blowing  air  was  provided  by  a  Carrier 
three-stage  centrifugal  compressor,  driven  by  a  300  horse- 
power General  Electric  induction  motor.   It  was  capable  of 
supplying  1900  cubic  feet  per  minute  at  29.5  psia.   Each 
plenum  was  supplied  and  metered  independently  by  a  gate  valve 
in  the  three-inch  supply  pipe.   Each  pipe  had  an  orifice 
plate  with  a  1.8  inch  hole  diameter  for  use  in  mass  flow 
measurements . 

The  flat  plate,  constructed  of  a  slab  of  phenolic  material 
24-inches  wide,  60-inches  long,  and  \   inch  thick  had  a  rounded 
leading  edge  in  order  to  avoid  leading  edge  separation.   The 
last  five  inches  was  tapered  and  hinged.   For  this  study  the 
hinged  portion  remained  at  zero  angle  of  incidence  with  re- 
spect to  the  flat  plate.   The  plate  was  mounted  at  the  center 
line  of  the  test  section  and  ran  from  a  point  7.5-inches  from 
the  leading  edge  of  the  plugs  to  a  point  five-inches  back  to 
the  end  of  the  plugs  as  shown  in  Fig.  III. 7. 

Figure  III.  9  shows  a  schematic  of  the  test  setup. 

C.   INSTRUMENTATION  AND  CALIBRATION 
1 .   Freestream  Sensors 

A  conventional  pitot  static  tube  located  22.5  inches 
aft  of  the  model  leading  edge  and  midway  between  the  flat 
plate  and  the  upper  model  section  was  used  to  measure  mean 

23  '  ' 


freestream  dynamic  pressure.   A  hot  wire  probe  2  9.5  inches 
aft  of  the  model  leading  edge  and  midway  between  the  plate 
and  the  lower  model  section  was  used  to  determine  magnitude 
of  velocity  perturbations  during  oscillations.   The  probe 
was  connected  to  a  locally  .manufactured ,  single  channel, 
hot  wire  anemometer  which  was  connected  to  a  Tektronix  555 
Dual-beam  oscilloscope  for  display  and  read  out.   Previous 
experimentation  [Ref.  12]  has  indicated  that  the  hot  wire 
circuitry  is  linear  with  freestream  velocity.   A  short 
calibration  run  was  made  which  verified  this. 

The  frequency  of  oscillation  was  obtained  from  a 
magnetic  pickup  mounted  out-board  of  the  uppermost  shutter 
blade  shaft,  as  seen  in  Fig.  III.  2.   The  output  frequency 
was  read  on  a  Dynascience  digital  counter. 

The  pressure  distribution  in  the  tunnel  was  measured 
by  a  series  of  flush  pressure  ports  one  inch  apart  in  the 
flat  plate.   They  ran  from  a  point  five  inches  in  front  of 
the  diverging  section  of  the  model  to  2.3  inches  from  the 
end  of  the  model.   The  ports  were  connected  to  a  40  tube, 
2  50  centimeter  upright  manometer  board.   The  median  reading 
of  the  oscillations  was  the  recorded  pressure  value. 
2 .   Blowing  Coefficient 

The  mass  flow  of  blowing  air  was  measured  by  a  pair 
of  orifice  plates  in  the  three  inch  supply  lines.   The  ori- 
fice plates  were  calibrated  in    situ.       This  was  accomplished 
by  connecting  the  blower  supply  lines  to  a  76.5  cubic  feet 
per  minute  rotameter  manuf actured  by  Fischer  and  Porter 
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Company.   The  pressure  drop  across  the  plates  was  measured 
by  a  36  inch  U-tube  water  manometer.   The  results  of  the 
calibration  were  plotted  against  the  ASME  values  as  shown 
in  Fig.  III. 10.   The  differences  were  so  small  as  to  be 
negligible  ,  therefore  it  was  deemed  reasonable  to  extrapolate 
the  curve  along  the  same  form  as  the  ASME  curve  to  cover 
values  not  covered  by  the  rotameter.   Figure  III. 10  was  then 
used  to  find  the  volume  flow  of  air  for  any  value  of  pres- 
sure drop. 

The  pressure  downstream  of  the  orifice  plate  was 
measured  with  a  36  inch  U-tube  mercurcy  manometer  vented  to 
the  atmosphere.   Air  temperature  remained  essentially  con- 
stant from  the  compressor  to  the  plenums.   With  the  volume 
flow,  pressure  and  temperature,  the  mass  flow  could  easily 
be  calculated. 

The  blowing  jet  velocity  for  calculation  of  C   has, 
historically,  been  defined  as  the  velocity  at  the  exit  as- 
suming isentropic  from  some  pressure  and  temperature  inside 
the  supply  duct.   These  were  measured  with  a  special  purpose 
Kiel  temperature  probe.   The  pressure  output  was  connected 
to  the  same  upright  manometer  board  as  the  pressure  distri- 
bution ports.   The  thermocouple  for  temperature  measurement 
was  connected  to  a  Leeds  and  Northrup  portable  potentiometer. 
3 .   Visualization  of  Separation 

Separation  lines  were  visualized  with  soft  twine 
tufts  placed  along  the  diverging  planes.   The  lines  of  tufts 
spanned  the  test  section  to  within  three  inches  of  each 
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sidewall  and  were  spaced  at  one  inch  intervals  across  the 
test  section  and  two  inch  intervals  down  the  plane  as  seen 
in  Fie.  III. 11. 
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Figure  III. 3.   Photograph  of  Shutter  Valve  Drive  Mechanism. 
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Figure  III. 6.   Photograph  of  Test  Section  and  Control  C 
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Figure  III. 11.   Photograph  of  Diverging  Planes  from  Downstream. 
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IV.   TEST  PROGRAM 

A.   EXPERIMENTAL  PROCEDURE 

The  basic  test  program  was  to  find,  for  some  value  of 
freestream  dynamic  pressure ,  the  value  of  blowing  coefficient 
that  was  just  sufficient  to  maintain  an  attached  boundary 
layer,  then  to  study  the  effects  of  varying  the  frequency 
of  oscillation  on  the  blowing  coefficient  required. 

The  first  problem  was  to  select  values  of  mean  free- 
stream  dynamic  pressure  for  the  study.   The  dynamic  pres- 
sures chosen  were  those  which  indicated  5,  7.5,  and  10  cen- 
timeters on  the  micro-manometer.   This  decision  was  based 
on  three  criteria: 

1.  Below  5  centimeters  the  atmospheric  wind  conditions 
at  the  inlet  began  to  affect  the  test  section  conditions. 

2 .  Above  10  centimeters  the  fans  could  not  maintain  the 
pressure  at  higher  oscillation  frequencies.   This  could  be 
remedied  by  opening  the  fan  housing  and  altering  the  vane 
direction  but  this  would  have  vastly  increased  the  time  re- 
quired for  runs  and  would  have  had  an  effect  on  test  section 
conditions . 

3.  These  three  settings  were  easily  read  and  maintained 
during  tunnel  operations. 

These  three  dynamic  pressures  gave  values  of  mean  freestream 
velocity  of  about  96,  118,  and  138  feet  per  second,  respec- 
tively, although  this  varied  slightly  with  oscillation 
frequency . 
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The  initial  tests  were  dictated  by  the  wind  tunnel  set- 
up from  previous  work.   The  set  of  belts  and  pulleys  con- 
necting the  electric  drive  motor  to  the  shutter  valve,  at 
the  beginning,  corresponded  to  a  frequency  range  of  2.4  8  to 
14.4  cycles  per  second.   Then  on  the  basis  of  the  previous 
test  results  the  range  was  expanded  outward  to  cover  fre- 
quencies of  1.6  to  16  cycles  per  second.   (See  Results  for 
detailed  explanation.) 

The  detailed  test  procedure  was  as  follows : 

1.  Start  tunnel  and  stabilize  at  desired  dynamic 
pressure . 

2.  Adjust  blowing  to  barely  maintain  attached  flow 
down  the  diverging  plane. 

3.  Take  data  for  pressure  distribution,  blower  plenum 
conditions  and  mass  flow  of  blown  air. 

4.  Start  shutter  valve  and  set  to  the  desired  frequency, 
Adjust  dynamic  pressure  as  necessary. 

5 .  Repeat  step  2 . 

6.  Take  data  for  magnitude  of  oscillation. 

7.  Repeat  step  3. 

8.  Adjust  shutter  valve  for  new  frequency  and  repeat 
steps  5,  6,  and  7  in  order. 

B.   DATA  REDUCTION 

The  data  reduction  was  accomplished  using  a  Fortran 
language  computer  program  on  the  IBM  360/67  computer  located 
at  the  Naval  Postgraduate  School.   The  program  is  shown  in 
a  separate  section  beginning  on  page  61. 
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The  calculations  involved  are  of  mean  free stream  static 
pressure,  density,  velocity  and  magnitude  of  the  oscillation 
of  mean  freestream  velocity  in  percent.   Then  the  pressure 
distribution  was  non-dimensionalized  with  respect  to  free- 
stream  dynamic  pressure.   Finally  the  blowing  coefficient 
was  calculated  for  the  upper  and  lower  blowers. 

In  addition  the  computer  program  plotted  the  pressure 
distribution  for  constant  frequency  and  blowing  coefficient 
versus  frequency  for  constant  freestream  dynamic  pressure. 
The  pressure  distribution  plots  were  used  for  further,  manual, 
data  reduction.   The  blowing  coefficient  variation  with  fre- 
quency was  a  guide  to  further  experimentation  in  the  early 
stages,  and  for  presentation  of  results. 
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V.   RESULTS 

The  results  of  this  study  are  presented  in  Figures  V.2 
through  V.ll. 

One  of  the  greatest  difficulties  encountered  during  ex- 
perimentation was  establishing  a  reasonable  definition  of 
separation  in  an  oscillating  flow.   Turbulent  oscillating 
flow  does  not  lend  itself  to  a  single,  unambiguous  defini- 
tion of  separation  because  the  flow  tends  to  separate  pro- 
gressively.  Flow  areas  may  be  identified  that  are  attached 
always ,  separated  always ,  and  alternately  attached  and 
separated  during  a  single  cycle  of  oscillation.   The  method 
of  visualizing  separation  with  rows  of  tufts  led  to  the 
definition  used  in  this  study.   It  was  assumed  that  separa- 
tion had  occurred  if  the  flow  reversed  itself  within  the 
boundary  layer  during  any  part  of  a  cycle.   The  definition 
used  in  this  study  does  not  coincide  with  Despard's  defini- 
tion [Ref.  12]  for  separation  in  laminar  boundary  layers. 
Despard  proposed  separation  as  commencing  with  the  initial 
occurrence  of  zero  velocity  or  reverse  flow  throughout  the 
entire  cycle.   Despard's  definition  cannot  be  easily  used 
for  flow  visualization  studies.   Preliminary  testing  re- 
vealed that  below  about  1.6  cycles  per  second  the  flow  acted 
in  a  quasi-steady  manner.   That  is,  the  blowing  required  to 
maintain  attachment  was  dependent  on  the  instantaneous 
velocity.   Therefore  it  was  necessary  to  provide  blowing 
equal  to  that  required  by  the  maximum  magnitude  of  the  adverse 

'il 


pressure  gradient  during  a  cycle.   Also,  above  16  cycles 
per  second,  the  oscillations  were  so  rapid  that  it  became 
impossible  to  detect,  by  eye,  the  onset  of  separation  with 
the  method  in  use.   In  fact,  as  the  frequency  approached  16 
cycles  per  second  the  decision  as  to  whether  to  flow  was  or 
was  not  separated  became  more  arbitrary  because  the  tufts 
would  oscillate  at  a  frequency  close  to  tne  limit  of  the 
eyes  ability  to  discern  change.   For  the  reasons  cited  above 
the  study  was  continued  to  the  frequency  range  from  1.6  to 
16  cycles  per  second. 

One  other  problem  encountered  with  the  experimental  setup 
was  that  of  three  dimensional  flow  in  the  test  section.   It 
has  been  suggested  that  the  study  of  separation  on  a  flat 
plate  in  a  turbulent  oscillating  flow  with  strong  adverse 
pressure  gradient  could  be  accomplished  by  using  a  model 
similar  to  that  used  in  this  study.   To  achieve  this  it 
would  be  necessary  to  have  the  flow  attached  across  the 
entire  plane  of  the  diverging  section.   This  should  allow 
the  boundary  layer  growing  on  the  flat  plate  to  become  de- 
tached due  to  the  impressed  pressure  gradient.   Tufts  were 
placed ,  therefore ,  on  both  the  upper  and  lower  surfaces  of 
the  flat  plate,  as  well  as  along  the  diverging  plane.   Boun- 
dary layer  growth  in  the  corners  and  on  the  sidewalls  caused 
a  wedge  shaped  area  of  separation  to  appear  on  the  diverging 
plane  as  shown  in  Figure  V.l.   Tufts  close  to  the  sidewalls 
indicated  this  corner-sidewall  separation  but  those  close  to 
the  centerline  of  the  tunnel  in:!'   ted  the  flow  was  still 
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attached  to  the  diverging  plane.   It  was  necessary,  therefore, 
to  observe  only  tufts  in  the  center  of  the  test  section  dur- 
ing this  study.   The  flow  was  considered  attached  if  at  least 
five  contiguous  tufts  in  one  row  indicated  that  the  flow  was 
attached.   This  definition  was  somewhat  arbitrary,  but  its 
consistent  application  led  to  some  meaningful  trends.   Some 
separation  did  occur  from  the  bottom  side  of  the  plate  when 
the  mean  freestream  velocity  was  138  feet  per  second  with 
high  blowing  coefficient.   Unfortunately,  this  separation 
was  not  reproducible.   The  flow  never  did  detach  from  the 
top  side  of  the  plate.   These  observations  show  that  it  will 
be  necessary  to  blow  along  the  sidewalls  to  prevent  the  cor- 
ner separation,  thereby  ensuring  more  predictable  separation 
characteristics  on  the  plate. 

To  analyze  the  results  it  was  necessary  to  determine  the 

parameters  on  which  C    (the  minimum  blowing  required  to 

yR 
maintain  attached  flow)  is  dependent.   It  was  assumed  that 

the  jet  velocity  required  V.   to  attach  the  flow  was  a  func- 

DR 
tion  of  the  density,  p;  mean  freestream  flow  velocity,  U^ ; 

the  length  of  the  plane  on  which  attachment  is  measured,  £; 
viscosity,  u;  the  magnitude  of  the  perturbation  of  the  free- 
stream  velocity,  AU  ;  the  frequency  of  the  oscillation,  co ; 
the  slot  height,  h;  and  the  effect  of  the  pressure  gradient 
as  measured  by  APF«6,  where  6  is  the  boundary  layer  thickness 
at  the  diffuser  exit  and  APp  the  pressure  difference  from 
the  reference  or  undisturbed  condition  to  the  exit.   This 
dependence  is  expressed  as  Equation  V-l. 
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v^   =£,EP»U  *,AP  .6,p,AU  u.h].  (V-l) 

JR 

Bernoulli's  equation  may  be  used  to  express  the  pressure 

coefficient  as 

APE  UE 

^^  oo  oo 

where  Up  is  the  velocity  at  the  exit  of  the  diverging  section. 
Substituting  V-2  into  V-l  and  non-dimensionalizing  with  re- 
spect to  p,  U^,  and  I    leads  to 

v  (V-3) 

iD  -  Ul  AU 

(  U—  )  =  f^C  J-         )-[l  -         ^         ],(  ^         ),(  jp         ),(  jp-         ),(  J  )} 

OO  CO  r       OO  CO  oo 

It  is  convenient  to  replace  u/pU  £  by  its  reciprocal,  the 
Reynolds  number. 

C.   may  then  be  expressed  as 


yR 


m.V.     pbhV? 

3  JR       ]R 


'uR    JspU2  3    ispU2(b£) 


(V-4) 


where  b  is  the  width  of  the  slot.   This  may  be  simplified  to 

V. 

C    =  2(  y  )(  ~-    )2.  (V-5) 

yR       *      U» 

From  equation  V-3  and  V-5  the  critical  blowing  coefficient 

may  be  written  as 

C    =  f3af).Cl-  gf  ],(  —    >,(£*),(  —  ),(  £  )>(V-6) 

K  J^  oo  r  CO  OO 


Equation  V-6  shows  that  C    is  a  function  of  five  independent 

yR 
parameters.   In  the  experimental  study  the  pressure  gradient 

parameter 

i+4 


«,  a  ;u   u2  J 

00 

was  determined  by  the  model  and  tunnel  configuration.   Al- 
though the  pressure  gradient  was  measured  (see  later  dis- 
cussion) there  was  no  measurement  or  control  of  6  possible 
in  this  study.   This  parameter  may  be  most  important  and 

further  studies  should  be  made  while  observing  6 .   Three 

5  5 

Reynolds  numbers  were  considered;  5.9  3x10  ,  7.31x10   and 

5 

8.5  3x10  .   These  Reynolds  numbers  corresponded  to  test  sec- 
tion freestream  velocities  of  96  f ps ,  118  fps,  and  138  feet 
per  second  (dynamic  pressures  of  5.0,  7.5,  and  10.0  centi- 
meters of  water),  respectively. 

Oscillation  frequency  was  varied  continuously  between 

1.6  and  16  cycles  per  second  which  resulted  in  a  reduced 
frequency  range  of  approximately  .1  to  1.1.   The  percent 
perturbation  of  velocity  AU  /U   was  a  function  of  frequency, 
shutter  valve  plate  size,  and  freestream  dynamic  pressure. 
It  was  not  independently  controllable  in  this  study  since  a 
single  plate  size  was  used.   AU  /U  ,  therefore  was  dependent 

or"  oo    oo  *  *■ 

on  frequency  and  dynamic  pressure.   Since  AU^/U^  did  vary  as 

the  reduced  frequency  was  changed,  the  resulting  frequency 

variation  in  C    measured  could  not  be  determined  at  constant 

yR 
values  of  all  the  correlation  parameters. 

Reference  [18]  states  that  experiments  have  shown  that 

for  steady  flows  C    is  not  a  function  of  h/£  for  high  slot 

MR 
Reynolds  Numbers,  V.h/v.   Therefore,  h/l    is  normally  omitted 

from  Equation  V-6 .   Blowing  characteristics  could  vary  sig- 
nificantly with  slot  height  at  low  values  of  the  slot  Reynolds 

M5 


numbers.   Slot  Reynolds  numbers  encountered  in  this  study  varied 

from  6.06xl03  to  11.56xl03  for  the  upper  slot,  and  6.42xl03  to 

3 
11.74x10   for  the  lower  slot.   To  independently  check  the  re- 
sults of  Ref.  [18]  for  the  case  of  unsteady  flows,  this  study 

-3  -3 

used  two  non-dimensional  slot  heights  of  3.15x10    and  4.23x10 

Figures  V-2,  3  and  4  clearly  indicate  that  frequency  has 

a  definite  effect  on  blowing  requirements  to  maintain  an  at- 

5 
tached  flow.   At  a  Reynolds  number  of  5.93x10  ,  figure  V-2 


shows  that  there  is  no  apparent  change  in  C    from  the  non- 

yR 
oscillating  condition  to  a  reduced  frequency  of  0.5  where 

C    increases  by  approximately  30  percent.   The  required 

yR 
flowing  gradually  decreases  through  the  remainder  of  the 


frequencies  tested  in  this  study.   For  the  Reynolds  number 

of  7.33x10  -  Figure  V-3.  the  initial  values  of  C    at  a 

reduced  frequency  of  0.1  are  approximately  40  percent 

above  the  non-oscillating  condition.   At  a  reduced  frequency 

of  0.13,  C    drops  to  approximately  30  percent  greater.   At 

yR 
0.48  the  blowing  required  increases  to  as  much  as  50  percent 

above  the  non-oscillating  condition  after  which  it  drops 
back  down  to  approximately  30  percent  for  the  higher  fre- 
quency oscillations.   Figure  V-4  shows  the  C    behavior  for 

5  Pr 

a  Reynolds  number  of  8.53x10  .   It  shows  that  a  reduced 

frequency  of  0.08,  C    is  approximately  40  percent  greater 

yR 
than  the  non-oscillating  case.   C    rises  rapidly  to  100 

yR 
percent  higher  at  a  reduced  frequency  of  0.12.   C    stays  at 

yR 
this  level  until  a  reduced  frequency  of  0.29  is  reached 

where  it  drops  to  50  percent  above  the  non-oscillating 
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condition.   At  a  reduced  frequency  of  0.4,  the  blowing  re- 
quired then  rises  rapidly  to  70  percent  above  the  non-oscil- 
lating condition  followed  by  a  gradual  decrease  through  the 
remainder  of  the  frequencies  studied. 

There  was  a  significant  difference  between  C    for  the 

yR 
upper  and  lower  blowers.   This  was  not  expected,  in  spite 

of  the  differences  in  slot  heights,  as  indicated  in  pre- 
vious discussion,  because  C   is  a  momentum  coefficient. 
Therefore,  although  the  upper  blower  had  less  mass  flow, 
its  greater  jet  velocity  should  produce  the  same  attachment 

effectiveness  as  the  lower  blower.   The  upper  blower  C 

yR 
was  approximately  15  percent  higher  than  that  of  the  lower 

5  5 

blower  at  R-  of  5.93x10  .   For  R„  of  8.53x10   the  difference 

5 
was  20  percent.   But  at  R{  equal  to  7.31x10  ,  the  upper 

blower  required  approximately  4  percent  less  blowing  than 

the  lower  blower.   This  inconsistency,  with  test  section 

velocity,  of  the  difference  between  upper  and  lower  C 

yR 
might  indicate  that  the  upper  and  lower  diverging  sections 

experience  different  flow  environments .   Since  separation 

is  intimately  linked  with  both  previous  boundary  layer 

growth  and  disturbing  influences ,  the  C    differences  may 

yR 
be  explained  by  the  presence  of  some  downstream  flow  asym- 
metry that  is  a  function  of  test  section  velocity. 

It  was  recognized  that  the  variation  in  C    was  depen- 

yR 
dent  on  the  pressure  gradient  and  the  magnitude  of  the 

perturbation  of  freestream  velocity,  both  of  which  exhibit 

variations  with  frequency  in  this  test  setup.   The  median 
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pressure  coefficient,  C  ,  on  the  flat  plate  was  plotted  versus 

position,  x,  along  the  plate.   x/£  equal  to  zero  corresponds 

to  the  blowing  slot  location.   Typical  pressure  distributions 

are  shown  in  Figure  V-5  for  the  cases  of  steady  flow  with  no 

blowing,  steady  flow  with  blowing,  and  oscillating  flow  with 

blowing  respectively.   The  pressure  gradient  was  determined 

by  estimating  the  slope  of  the  best  straight  line  fit  of  the 

plotted  pressures.   The  resulting  slope  was  plotted  against 

the  reduced  frequency  for  each  of  the  freestream  velocities 

run.   These  data  showed  similar  variation  of  C    versus  re- 

Px 
duced  frequency  and  so  are  plotted  together  in  Figure  V-6 . 

An  examination  of  Figure  V-6  shows  that  the  pressure  gradient 

variation  with  frequency  was  not  the  primary  cause  of  the 

exhibited  large  variations  of  C    with  f requenuy .   Similar 

UR 
plots  of  percent  perturbation  in  freestream  velocity,  as 


measured  by  the  hot  wire,  versus  reduced  frequency  are  shown 

5 
in  Figures  V-7,  8  and  9.   The  cases  of  R„  equal  to  5.93x10 

and  7.31x10   show  a  gradual,  almost  linear,  decrease  of 

perturbation  velocity  with  frequency  but  the  R„  of  8.5  3x10 

case  showed  the  characteristics  seen  in  Figure  V-9.   There 

was,  however,  no  indication  that  these  variations  in  the 

velocity  perturbation  were  responsible  for  the  shape  of  the 

large  variation  of  C    with  frequency  seen  in  Figure  V-2 ,  3 

VlR 
and  4. 

The  blowing  requirement  to  maintain  attached  flow  has 
been  shown  to  be  frequency  dependent.   This  dependence  ex- 
hibits characteristics  which  suggest  resonant  behavior. 
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This  resonance  might  be  the  result  of  flow  confinement. 
Wave  reflections,  due  to  the  tunnpl  walls,  could  affect  the 
conditions  in  the  test  section.   These  reflected  waves  could 
either  reinforce  or  dissipate  the  shutter  valve  induced 
pressure  pulses.   Another  source  of  resonance  lies  in  the 
physical  makeup  of  the  tunnel.   It  was  observed  during  the 
experimental  runs  that  certain  frequency  bands  at  each  tun- 
nel test  section  velocity  made  the  wind  tunnel  vibrate  a 
great  deal,  sometimes  causing  the  test  section  to  be  travel- 
ing as  much  as  two  inches  longitudinally.   These  resonant 
frequency  bands  are  shown  in  Figures  V-2,  3  and  4.   This 
resonant  behavior  of  the  tunnel  occurred  around  a  frequency 
of  approximately  8.7  cps.   Tunnel  vibrational  modes  are  seen 
to  be  influencing  the  blowing  characteristics  observed.   Dif- 
ferent test  configurations  may  exhibit  different  characteristics 
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Figure  V.l, 


Photograph  of  Tuft  Showing  Corner  Sidewall 
Separation. 
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Figure  V.3.   Blowing  Coefficient  Required  vs.  Frequency, 
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Figure  V.9.   Velocity  Perturbation  vs.  Frequency,  R   =  8.5  3x10 
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VI.   CONCLUSIONS 

From  the  results  obtained,  the  following  conclusions  may- 
be drawn: 

1.  There  is  a  definite  effect  on  the  blowing  required 
to  maintain  an  attached  boundary  layer  in  a  strong  adverse 
pressure  gradient  as  oscillations  are  superimposed  on  the 
freestream  velocity  and  the  frequency  is  varied. 

2 .  The  frequency  dependence  of  the  blowing  requirements 
exhibits  characteristics  which  suggest  resonant  behavior. 

3.  Considerable  three-dimensional  flow  is  produced  in 
the  present  test  setup  and  sidewall  blowing  will  be  required 
to  produce  two-dimensional  separation  on  the  plate. 
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VII.   REC  OMMHttDAJg  ON  S 

1.  C    is  a  function  of  the  five  parameters  discussed 
in  the  Results  section.   More  tests  should  be  made  in  which 
only  one  of  these  is  varied  at  a  time.   As  discussed  earlier, 
measurement  of  the  boundary  layer  thickness  is  most  desirable 

2.  The  velocity  profiles  throughout  the  flow  field  need 
to  be  studied  as  a  function  of  time  while  varying  oscillation 
frequency.   To  accomplish  this  a  method  for  flow  measurement 
must  be  developed  which  will  fulfill  the  following  require- 
ments;  non-interference  with  the  flow,  capable  of  measuring 
turbulent  variations ,  capable  of  scanning  a  relatively  large 
section  of  the  flow,  and  capable  of  resolving  rapid  changes 
in  flow  velocities  at  high  oscillation  frequencies. 

3 .  An  investigation  of  the  natural  frequencies  of  the 
tunnel-model-air  flow  configuration  should  be  made  to  deter- 
mine the  effects  on  pressure  gradient,  magnitude  of  oscil- 
lating velocity  perturbation,  and  blowing  requirements  to 
maintain  an  attached  boundary  layer. 

4 .  Sidewall  blowing  should  be  used  to  ensure  two  dimen- 
sional separation  of  the  flat  plate  for  further  studies. 
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CCMPUTER    PROGRAM 


Q  ********************************************************* 

c  *                                                                           * 

r  if                                                                                                                                                                                  # 

C  *       THIS    PROGRAM    CALCULATES    THE    FOLLOWING     .                                        * 

C  *               1.       MEAN    FREESTREAM    VELOCITY    IN    FEET     PER    SECOND         * 

C  *               2.       OSCILLATION    IN    RADIANS    PEP    SECOND                                     * 

C  *               3.        FREESTREAM    VELOCITY    PERTLRBITION     IN    PERCENT         * 

C  *               h,       DIMENSIONLESS    PRESSURE    DISTRIBUTION                               * 

C  *                          A.       ALSO    PLOTTED    ON    CALCCMP    PLOTTER    WITH    THE       * 

C  *                                     ABSCISSA    THE    X    POSITION    IN    INCHES                          * 

C  *                                     BEGINNING    AT    THE    UPSTREAM    POSITION    AND            * 

C  *                                     THE    ORDINATE    IS    THE    PRESSURE    COEFFICIENT      * 

C  *               5.       JET    VELOCITY    OF    EJECTED    AIR    FOR    BOTH    UPPER            * 

C  *                          AND    LOWER    BLOWERS    IN    FEET    PER    SECOND                            * 

C  *               6.        MASS    FLOW    OF    EJECTED    AIR    IN                                                     * 

C  *                          POUNDS    MASS/SECOND                                                                              * 

C  *               7.        BLOWING    COEFFICIENT    FOR    UPPER    AND                                     * 

C  *                          LOWER    BLOWERS                                                                                            * 

C  *                          A.       PLOTTED    ON    CALCOMP    PLCTTER    WITH    THE                    * 

C  *                                     ABSCISSA    THE    FREQUENCY    IN    RADIANS    AND               * 

C  *                                     THE    ORDINATE    THE    BLOWING    COEFFICIENT                  * 

C  *       OUTPUT     IS    LINE    PRINTED,     SEPARATED    BY    FREQUENCY                    * 

C  *       AND    IS    LABELLED                                                                                                         * 

c  *                                                                                                                                     ^ 

C  *                                                                                                                                                 * 

C.'.  .',  O.  %'.   »'.  v',  O,  .1.   ,'-  .'.  ,'.  J.  ,!.  >'  ■  J,  ,'.  .',  »'    .'.-  ,''  ■'-  »<*  v'-  »'-  V.  n'-  ••-  »'/  J-  -I.  .'   .1-  .'.  v1.  *'■  ^',  >'*  J-  ^  4-  J,  J,  *'.  Um    . .".  -J'     -'-  -''  -  -  -'*    -"-  -->-   .'.-■-  »'-  %>„  .'. 
*•»  *i-  *•»  -I*  -*>  A*  «v  *V»  3f»  *»■•  **  1"  A*  'f  ~r-  *V  *Y"  *V*  *»»*»-■(»  <V   -I*  »•*  f  *(*  •*<*  *V  f  ">*  "**  A'  't-   *Y*  *+■  *\  •  'f  'f   -I*  *t~   T"  T>  T1  "t*  T*  *"*  *r  *1*  "**  T  Ti*  ***  'f  1*  TT  If  -T- 
c 
c 
c 

C  INPUT  -REQUIRED 

C  PATM    =    ATMOSPHERIC    PRESSURE     (INCHES    OF    MERCURY) 

C  T    =    OUTSIDE    AIR    TEMPERATURE     (DEGREES    FAREHEIT) 

C  FREQ    =     FREQUENCY    (CYCLES    PER    SECOND) 

C  Q    =    MEAN    FREESTREAN    DYNAMIC    PRESSURE     (CENTIMETERS 

C  OF    WATER) 

C  QMAX    =    MAXIMUM    MAGNITUDE    OF    VELOCITY    OSCILATIONS 

C  (VOLTS    OUTPUT    CF    HOTWIRE    ANEMOMETER) 

C  QMIN    =    MINIMUM    MAGNITUDE    OF    VELOCITY    OSCILATIONS 

C  (AS    ABOVE) 

C  PAMB    =    AMBIENT    PRESSURE     (CENTIMETERS    CF    WATER, 

C  REFERENCE    HEIGHT    FROM    MULTI-TUBE    MANOMETER    BOARD 

C  P    =    PRESSURE    DISTRIBUTION,     16    VALUES    PUT     IN    1X16 

C  VECTOR    (CENTIMETERS    OF    WATER) 

C  P2U    =    PRESSURE    JUST    DOWNSTREAM    OF    THE    ORIFICE    PLATE 

C  .     IN    THE    SUPPLY    LINE    TO    THE    UPPER    BLOWER 

C  (CENTIMETERS    OF    WATER) 

C  P2L    =    PRESSURE    AS    ABOVE     BUT    LOWER    BLOWER 

C  VFLOU    =    VOLUME    FLOW    OF    PLOWING    AIR    TO    UPPER    BLOWER 

C  (CUBIC    FEET    PER    MINUTE) 

C  VFLOL    =    VOLUME. FLOW    AS    ABOVE    FOR    LOWER    BLOWER 

C  TPU    =    PLENUM    AIR    TEMPERATURE,     UPPER    BLOWER    PLENUM 

C  (DEGREES    FARENHEIT) 

C  TPL    =    PLENUM    AIR    TEMPERATURE,     LOWER    BLOWER 

C  PPU    =    PLENUM    AIR    PRESSURE,     UPPER    BLOWER 

C  (CENTIMEPS    OF    WATER    DIFFERENCE    FkOM    ATMOSPHERIC) 

C  PPL    =    PLENUM    AIR    PRESSURE,     LOWER    BLOWER 

C 

C  INPUT    GOES    ON    FOUR    DATA    CARDS,     THE    FIRST    HAS    7    NUMBERS 

C  THEY    ARE,     PATM,     T,     FREQ,    Q,     QMAX,     GMIN,     PAMB. 

C  THE    NEXT    TWO    CARD    HAVE    THE     16    NUMBERS    OF    PRESSURE 

C  DISTKIBUTIuN.        THE    LAST    CARD    HAS,     P2U,     VFLOL,     TPU, 

C  PPU,     P2L,     VFLOL,     TPL,     PPL.       THE    NUMbtRS    ARE    INPUT 

C  AS    F10.4 

Gl 


c 

c 

C      LINE  OUTPUT 

C  FREQUENCY  (CYCLES  PER  SECOND) 

C  Q  (CENTIMETERS  OF  WATER) 

C  VELOCITY  (MEAN  FREFSTREAM  IN  FEET  PER  SECOND) 

C  PERCENT  (PERTURBATION  OF  VELOCITY) 

C  OMEGA  (OSCILATICNS  IN  RADIANS  PER  SECOND) 

C  CMU  (BLOWING  COEFFICIENT,  UPPER  BLOWER) 

C  VJU  (JET  VELOCITY,  UPPER  BLOWER  IN  FEET/SECOND 

C  FLOU  (MASS  FLOW  OF  JET,  UPPER  BLOWER  IN 

C  POUNDS  PEP  SECOND) 

C  CML  (BLOWING  COEFFICIENT,  LOWER) 

C  VJL  (JET  VELOCITY,  LOWER) 

C  FLOL  (JET  MASS  FLOW,  LOWER) 

C  STATIONS  (PRESSURE  DI  STR  I  BUT  I  GN«.TAP  POSITIONS 

C  WITH  1  AS  THE  DOWNSTREAM   POSITION) 

C  CP  (PRESSURE  COEFFICEiMT  AT  STATION) 

C 

C 

C      PLOTTED  OUTPUT 

C  PRESSURE  COEFFICIENT  VS  X  POSITION  (WITH  1  BEING 

C  UPSTREAM) 

C  BLOWING  COEFFICIENT  VS  OMEGA  (PLOT  OF  CML  AND  CML 

C  VERSUS  OMEGA  FOR  EACH  INPUT  DECK) 

C 

c 
c 

DIMENSION  CP(16,40),CMU(40),QMEGA(40),F(40)  , 
1PCSL(16),P(16) ,CML(40) ,LP0S(16)  ,  C PI ( 16 )  , CA ( 40 ) 
10  FORMAT  (7F10.4) 
15  FORMAT  (8F10.4) 

20  FORMAT  ( • 0 • , T3 , • Q ' , T17 , « VELOC I TY « , T30, ■ PERC ENT • , T44 , 
1' OMEGA1  ,T58, 'CMU'  , T72 ,  • VJU '  ,T8  4 , 'FLOU' , T98,  'CML' 
ltTllO, 'VJL* iT123i 'FLOL'  ) 
25  FORMAT  ( « 0 • , 10F 13 . 5 ) 

30  FORMAT  ( • 0 • , Tl 9 , • STAT  I  ON '  , T37 ,  • CP • / ) 
35  FORMAT  (•  «,I20,F26.8) 
40  FORMAT  PI') 
4  5  FORMAT  (6A8) 
50  FORMAT  COLAST  =  «,I10) 
55  FORMAT  COTHE  FREQUENCY  IS  '  ,F10.5) 
REALMS  LABELl/'CML      •/ 
REAL*8  LA8EL2/'CMU      •/ 
REAL*8  LABELS/'  '/ 

REAL*8  TITLEK12) 
READ  (5,45)  TI7LE1 
REALMS  TITLE2U2) 
READ  (5,45)  TITLE2 
DO  1000  J=l,28 

READ  (5,10)  PATM,T,FREQ,Q,QMAX,QMIN,PAMB 
READ  (5,15)  P 
DC  100  K=l,16 
LPOS(K) =K 
PCSL(K)=LPOS(K) 
100  CONTINUE 

READ  (5,15)  P2U,VFL0U,TPU,PPU,P2L,VFL0L,TPL,PPL 
C 

C      CONVERT  TEMPERATURES  TO  RANKINE 
C 

TEMP=T+459.688 
TPL=TPL+459.688 
TPL'=TPU+459.683 
QA(J)=Q 
F(J)---FREQ 
C 

C      CALCULATE  OMEGA,  MEAN  FREESTREAM  VELOCITY,  AND 
C      PERCENT  PERTURBATION 
C 

OMEGA( J  )=FRE0*6. 283185308 

PINF=(?(14)+P( 15)+P(16))/3.0 

RHGM (2.0-PAMB) -PINT) *PATM*  .  04  120696200/ PAM  E7TEMP 
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c 
c 
c 


c 

c 
c 


c 
c 

c 


VEL=2.0  20916944D0*SCRT(Q/RH0) 
PERCEN=C(QMAX-QMIN)/(QMAX+QMIM) )*100 

NCN  DIMENSICNALIZE  PRESSURE 

DC  200  1=1,16 
CP(I,JI=( PINF-Pl I ) )/Q 
20C  CONTINUE 

CALCULATE  BLOWING  PARAMETERS 

RH02U=(P2U+PATM)*.0412  069621D0/TPU 

FLOU=VFLOU*RH02U/60.0 

VJU=DSQRT(545  4.5  0  724D0*TPU*(1.0-< ( ( 2. 0*P A MB )-PI NF ) 
1/{(2.0*PAMB)-PPU) )**. 2862 24 12 5 6 DO) ) 

CMU(J) =(FL0U*VJU)/(Q*4. 3322714 L7D0J 

RH02L= (P2l+PATM)*.0412  0696  21DO/TPL 

FLGL=VFLQL*RH02L/60.0 

VJL=DSQRT(5454.50724D0*TPL*(1.0-< ( (2.0*PAMB )-PINF J 
l/(  (2.0*PAMB)-PPL) )**. 2362  24  125 6D0)  ) 

CML(J)=(FLGL*VJL)/(Q*4. 3322714 17D0) 

LINE  OUTPUT 


F(  J) 


VEL,PERCEN 
,VJL,FLOL 


300 
.000 


C 
C 


WRITE  (6,40) 
WRITE  (6,55) 
WRITE  (6,20) 
WRITE  (6,25)  Q, 
lVJU,FLOU,Ci-,L(  J) 
WRITE  (6,30) 
DO  300  1=1, 16 

WRITE  (6,35)  LPOS(  I), CPU, J) 
CONTINUE 
CONTINUE 

PLOTTING  ROUTINES 


OMEGA ( J),CMU( J) 


NFTS=23 

EXSC=15.0 

YSCL=.006 

MX  =  0 

MY=0 

CALL  QRAW(NPTS,0MEGA,CML,1, 1,LABEL1,TITLE1 
1 ,  EX  SC , YSC  L , 0 , 0 , MX , MY , 9 , 15 , 0 , LAS  T ) 

WRITE  (6,50)  LAST 

CALL    DRAW(NPTS,0MEGA,CMU,3,2,LABEL2,TITLE1 
1,EXSC,YSCL,0,0,MX,MY,9,15,0,LAST) 

WRITE     (6,50)     LAST 

DO    2000    J=l,28 

DO    1100    1=1,16 

K=17-I 

CPK  I)=CP(K,  J) 
1100    CONTINUE 

CALL    DRAW( 16,P0SL ,CP1 , 0 , 2 , L ABEL3 , T ITLE2 , 
12.0,0.0,0,0,1,0,8 ,8,0, LAST) 

WRITE     (6,50)     LAST 
2000    CONTINUE 

STOP 

END 
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